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Abstract 

The vacuum neutrino mixing is known to exhibit an approximate ^ — t symmetry, which was 
shown to be preserved for neutrino propagating in matter. This symmetry reduces the neutrino 
transition probabilities to very simple forms when expressed in a rephasing invariant parametriza- 
tion introduced earlier. Applications to long baseline experiments are discussed. 
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I. INTRODUCTION 



The tremendous progress in the last decade has made it possible to pin down, with im- 
pressive accuracy, many of the fundamental parameters in the neutrino sector. A complete 
picture, however, is still not available. Chief among the missing information is the deter- 
mination of the Vi3 element of the neutrino mixing matrix V, which, in turn, is crucial in 
ascertaining the CP violation effects in the leptonic sector. Given that direct CP violations 
in the quark sector [l[ have been well-established and accurately measured, it is impera- 
tive, from both the theoretical and experimental points of view, to assess the corresponding 
situation in the leptonic sector. Another unsolved puzzle concerns the neutrino mass spec- 
trum, in that there are the possibilities of either the "normal" or "inverted" orderings. It is 
certainly important to settle this question. 

While the fundamental parameters refer to those in vacuum, it has been well-established 
(see, e.g., Ref.(2l-[l5l|) that they are modified when neutrinos propagate through matter, 
by giving the neutrino an induced mass, which is proportional to its energy and to the 
medium density. Indeed, in the analyses of the solar neutrinos, certain features of the data, 
such as the modification of the energy spectra from the original, can only be understood 
by the inclusion of matter effects. With the advent of long baseline experiments (LBL, 
for an incomplete list, see, e.g., Ref.(l6-23l|). the induced mass can actually be "tuned" 
by changing the neutrino energy {E). This provides a powerful tool which can be used to 
extract fundamental neutrino parameters from measurements. 

In this work, we will use a rephasing invariant parametrization which enables us to obtain 
simple formulas for the transition probabilities of neutrinos propagating through matter of 
constant density. It was shown earlier that these parameters obey evolution equations as a 
function of the induced mass. In addition, these equations preserve the approximate /i — r 
symmetry jl^, 2^ which characterizes the neutrino mixing in vacuum. Incorporation of 



the /i — r symmetry for all induced mass values results in a set of very simple transition 
probabilities P(z/a J-'p)- In general, these formulas offer quick estimates of the various 
oscillation probabilities, using the known solutions obtained earlier. As an example, we will 
analyze P{ve i^ti) in detail, emphasizing its dependence on the neutrino parameters. 



II. THE REPHRASING INVARIANT PARAMETRIZATION 



Neutrino oscillations, being lepton-number conserving, are described in terms of a mixing 
matrix whose possible Majorana phases are not observable. Thus it behaves just like the 
CKM matrix under rephasing transformations, which leave physical observables invariant 
"261 1. To date, however, such observables are often given in terms of parameters which are 



not individually invariant. So it seems that the use of manifestly invariant parameters may 



be more physically relevant. Two such sets are known to be \Vij\ [27|, |28| and VaiVj3jV*jV^i 



1 i „ „ I -J I r - a r-j ulj 

29|. Recently, by imposing the condition detV = +1 (without loss of generality), another 



set was found, given by [20, 30l^32 



ijk 



ViiV2jVzk — R 



ijk 



iJ, 



(1) 



where the common imaginary part can be identified with the Jarlskog invariant J |29| . Their 
real parts are labeled as 



(i?123, -R23I, -R312; -R132, -R213, -R321) = {Xi,X2, X3; Hi, 1/2, Us)- 



(2) 
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The variables are bounded by —1 < (xj, yj) < +1 with yj < Xi for any and satisfy two 
constraints: 

detV = (xi +X2 + X3) - (yi + 1/2 + y-i) = 1, (3) 

{xiX2 + X2X3 + X3X1) - {yiy2 + 2/2^/3 + ym) = 0. (4) 
Eq. (jlj), together with the relation 

J2 = X1X2X3 - ymys, (5) 

follow 2^ from (the imaginary and real parts of) the identity ri23r23ir3i2 = ri32r2i3r32i. 
Thus, flavor mixing is specified by the set {x, y) plus a sign, according to J = ±V^. This 
sign arises since the transformation V —^V* , corresponding to a CP conjugation, leaves the 
real part {x,y) of Tij^ invariant, but changes the sign of its imaginary part (J). Note that, 
using l^jP, a complete parametrization also requires four elements plus a sign. 
The parameters (x, are related to the rephasing invariant elements \Vij\^ by 



W 




(6) 



One can readily obtain the parameters (x, y) from W by computing its cofactors, which form 
the matrix w with uFW = (detl^)J, and is given by 

xi +yi X2 + y2 X3 + ys 
w = I X3 + y2 xi+y3 X2 + yi I . (7) 



X2 + y3 X3 + yi xi + y2 



The relations between (x, y) and 



utf ^ v^,Vp,v:y* (8) 



are given by (using K^iVaj - VajVpi = Y.-yk (^afi-yeijkV*,,): 

^aP^^ijk ^ai '^j3j ^•yk 

-yk 

= \^aj\'^\^l^i\'^ + ^(^al3-y(^ijkV*jV^iV*i.. (9) 

■yk 

The second term in either expression is one of the F's (r*'s) defined in Eq. ([T]). Also, by 
using the constraint in Eq. Re(n^^) can be expressed in terms of quadratics in {x,y), a 
result which will be used later in Tables I and II. 



III. EVOLUTION EQUATIONS AND THE /x - r SYMMETRY 



For neutrinos in matter (of constant density), it was shown [3l|, |32j that, as a function 
of the induced mass A = 2\/2GFneE, the neutrino parameters satisfy a set of evolution 
equations which are greatly simplified by using the (x, y) variables. It was found that 



dA 



Xi 



1,2,3) 



(10) 
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where Di are the eigenvalues of the Hamiltonian. Also, the evolution equations for all (xj, Uj) 
can be obtained and are collected in Table I of Ref . 3l|, |32| . Of particular interest for our 
purposes are the equations: 



din J 

dA 



■{xi - yi) + {x2 - y2) , -(a:2 -1/2) + (a^s -Z/s) 



+ 



D^-D2 
{xi - yi) - (x3 - y-i) 



Do - D. 



(11) 



and 





-yi) = 


X2 - 

Di- 


1/2 
D2 


X3 — 


1/3 


^^^"^^^ 


-^2) = 


Xi 

Di 


- yi 
-D2 


D2 


- 2/3 

-D3 


^^^"^^^ 


-1/3) = 


X2 
D2 


- 2/2 

-D3 




- yi 



(12) 



Note that the quantities Di — Dj and Xi — yi form a closed system under the evolution 
equations, independent of other possible combinations of these variables. 

There remain two more independent evolution equations, which may be chosen as those 
for {xi + yi). We define 

-^i Xi 



Then 



dfli 



E 



1 



f2j Xi ~\~ yi 



[6ij{QiXk — flk^i) — Sik{^iXj — QjXi 



(13) 
(14) 



It follows that 



— ^ijkii^iXj — QjXi) — {QiXk — QkXi))]. 

d 



dA 



{xi + ?/i) = 



(15) 
(16) 



if {xj + yj) = 0. This condition is equivalent to W2i = W^i, i.e., — r exchange sym- 
metry. Thus, the evolution equations preserve the fi — t symmetry, which was established 
(approximately) for neutrino mixing in vacuum. 

Another useful property of the evolution equations is to establish matter invariants. For 



instance 133^36 



d 

dA 



[X,X2XsAl,Al,Al,] = 



where Xi is defined in Eq. (IT^ and 



A,j=D,-Dj. 



(17) 



(18) 



(Also, d{JAi2A23A3i)/dA = 0, as mentioned before [3l|, |32|). In addition, there is a simple 
relation 

^ ^ E(X,-X,)A,, =1. (19) 

i>j 



2dA 
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Eqs. fll7p and fll9p are three-flavor generalizations of the two-flavor results 



32 



^(^Z/^') = 0, (20) 

±[{x + y)D] = -l, (21) 

where x = V11V22 = cos^ 0, y = V12V21 = — sin^ 9, D = m\ — mf, in the usual notation. 

The vacuum neutrino masses are known to be hierarchical, Sq/Aq ^ 1/32 ^ 1, (5o = 
7712— ml, '^0 = |'^3~''^2l- There are two possibilities, the normal hierarchy (mg S> ml ~ m"^), 
or the inverted hierarchy (mg <C ml ~ m"^). In matter of constant density, mf — )■ Z)j, 
which are A-dependent. For the case of normal hierarchy, there are two A- values where the 
levels "cross", at the lower resonance, A = Ai, [d{Di — D2)/dA]Ai = 0, and at the higher 
resonance, A = Ah, [d{D2 — D3)/dA]A^ = 0. From Eqs. f|T2l) . one finds that rapid variations 
occur only for A to be near Ai or Ah. Let us denote by {Aq, Ai, Ai, Ah, Ad) the values of 
A in vacuum [Aq = 0), at the lower resonance [Ai), in the intermediate range {Ai), at 
the higher resonance (Ah), and in dense medium (A^). Then, the solutions for {X,Y) are 
well-approximated [3l|, l32| by two-flavor resonance solutions. 



For < v4 < Ai, 

A21 = [p^A'-2qi6oA + 6',Y/\ 
Xi = ^[pi-{p^A-qi6o)/A2i], 

X2 = ^[pi + {p^A-qi6o)/A2i], 

= (X3)o, (22) 

where Ajj = Di — Dj in matter, Xi = Xi — yi, pi = {Xi + ^2)0, qi = {Xi — X2)o- Note that 
(Xi)o = 2/3, (^2)0 = 1/3, and XiX2Al^ = constant. 
For Ai<A< Ad, 

A32 = [plA'-2qhAd + A^Y/\ 
Xi = {Xi)i, 

X2 = ^[ph-{plA-qhA,)/A32], 

X3 = l[Ph+ {pIA - g,A,)/A32]. (23) 

Here, A = A - Ai, and ph = {X2 + ^3)^, qh = {X2 - ^3)^, A^ = (A32)i are taken at 
A = A > So. Note that (Xi), = 0, (X2). = 1, (X3). = (^3)0 = \Vis\l < 1- Also, X2X3Ai2 
is an invariant as A varies. Thus, the product X2X3 has a resonance behavior near A ~ A/j. 



Note also that the minimum of A32 is at (A32)mm — 2y IV13I0A0. 

To obtain A21 for Ai < A < Ad and A32 for < A < Ai, one first notes from Eq. ([10 
that dA2i/dA ~ X2 for high A. Thus, a direct integration leads to 



A 



1 



21 



A, + phA - {plA' - 2qhA,A + A^) 1/2] (24) 
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FIG. 1: The variation of Ajj for the normal hierarchy (left) and the inverted hierarchy (right). 
Both the numerical (solid lines) and the approximate (dashed lines) solutions are shown. The 
approximate analytical solutions are given by Eqs. (22-25). We adopt A = 0.02 and the current 
upper bound for jVisI, e ~ 0.17. 




for Ai < A < Ad, where Si 
low A gives 



(A2i)i. Similarly, a direct integration of dAs2/dA ~ X2 for 



A32 = Ao + - [5o - PiA - {p]A^ - 2qi5oA + 5lY'\ 



(25) 



The solutions for A31 in both regions of A are obtained from A31 = A32 + A21. Note that 
the solutions for {] < A < A^ and for Ai < A < A^, should agree for A = Ai. This condition 
leads to Aj ~ Aq — Ai and 5^ ~ Aj. 

For inverted hierarchy, the behaviors of Xj near Ai are given by the same Eq. (122|) . 
However, for A > Ai, there is no longer a resonance. Instead, all Xi change slowly, so that 
Xi ~ 0, X2 — 1, X3 ~ 0, for A > Ai. The solutions for u are obtained by A — > —A. Thus, 
there is a resonance behavior near Ah, for the inverted hierarchy scenario. Otherwise all the 
changes are small. 

The accuracy of the approximate formulas in Eqs. f l22tl25]) can be assessed by numerical 
integrations of the exact equations, Eqs. (fTU]) and (IT^ . To do that we write 



W 



,3 ^'^/ ,^ 3 



+ 27] 



6„ 



(26) 



where (e, 77, A,^) ^ 1 in vacuum, and W reduces to the tribimaximal [37| matrix when 
e = r/ = A = e = 0. 

It should be emphasized that the parameters (e, r], A,^) carry quite distinct behaviors as 
A varies, as shown in the following. Eqs. ([6]) and fl26|) give rise to 



^ = W23- W33 = {X2 + Z/2) - {xi + yi), 
and from 1^21 — W31, we have 

6A = 3(x3 + y-i) - 2{x2 + 2/2) - (a^i + yi)- 



(27) 



(28) 



With the constancy of Xj + yj, one concludes that ^ 
Wii + W12 = 1 - e^, we have 



A ~ as A varies. In addition, since 



dA 



_d_ 

"dA 



[{x2-y2) + (xi-yi)]. 



(29) 
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and 



— = 0, (for low A) 

— = -2{X2 - i/2)(a;3 - yi)l{D2 - D^), (for high A) (30) 



Furthermore, one obtains from Wu — Wu that 



^ = ^(1 - e') + \lix2 - y2) - ix^ - y,)], (31) 



and 



— = -(x2 - y2)(a;i - yi)/(£'i - A), (for low A) 
dn 2 

^ = 3(0:2 - 1/2)(X3 - y3)/{D2 - D,), (for high A) (32) 

Thus, 7] and can change considerably as functions of A, but A ^ ^ ~ throughout. 

For numerical integrations, Eqs. ([6]) and fl26l) suggest the following initial values in vac- 
uum: 

xio = ^(2 - 3A - 2e2), yio = ^(-2 - 3A + 26^), 
2:20 = ^(1 - 3A - e^), ^20 = ^(-1 - 3A + e^), 

D 

^30 = + 2/30 = ^(A-e2), (33) 

where ^ = 77 = is chosen and the terms in 0{Xe^) are ignored. We shall choose the initial 
values e = 0.17 and A = 0.02, which correspond to the experimental bounds |Vi3p < 0.03 
'3^ and an assumed CP violation phase cos ip = 1/4, respectively. The numerical solutions 



for the (x, y) parameters, the squared elements of the mixing matrix, and J in matter follow 
directly and are shown in Figs. 2-5 in Ref. j32|. Our choice of A 7^ signifies a small /i — r 
symmetry breaking, the solutions verify that {xi + yi) remain negligible for all A values. In 
addition, we show in Fig. 1 both the numerical and the approximate solutions for Aij in 
matter. Note that the hierarchical relation among the Aj^'s varies in matter and plays an 
important role in the oscillatory factor sin^ $jj of the probability functions. It is seen that 
A21/A31 ~ A21/A32 ~ 1/32 < 1 (normal hierarchy) and A21/A23 ^ A21/A13 ~ 1/32 < 1 
(inverted hierarchy) for < A ^ Aj. While in Ai < A ^ Ad, the Aj^'s are less hierarchical: 
A21/A31 ~ A21/A32 ^ 1/5 (normal) and A13/A21 ~ A13/A23 ~ 1/5 (inverted). 

IV. THE PROBABILITY FUNCTIONS 

The neutrino transition probability in matter is given by j38| 

P{^a ^ ^p) = 5^/3 - 4 5] M'n^f) sin^ 

j>i 

+ 2^Im(n:^/)sin2$,„ (34) 

j>i 
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complete 


with Xj + j/j — 






-X1X2 - X1X3 + xiy2 + yiys 


-2xiX2 


1 


^31 


X1X2 + x^yi - yiy2 - yws 


-2X1X3 


< 1 




-X1X2 - X2X3 + 2:2^3 + yiy2 


-2X2X3 


-2x2X3 


F2I 


+X1X3 + X2yi - yiy2 - yivs 


-2xiX2 


^ 1 


F31 


-X1X2 - X1X3 + xiy3 + yiy2 


-2X1X3 


< 1 


^32 


-X1X2 + X3y2 - yiy2 - ^22/3 


-2X2X3 


-2x2X3 




-X1X3 - X2X3 + X3y3 + yiy2 


-(X3/2) + X1X2 


-X3/2 




X1X3 + X2y2 - yiy2 - 2/22/3 


-(X2/2) + X1X3 


-X2/2 


F!i2 


-2:1X2 - X1X3 + Xiyi + 7/22/3 


-(xi/2) + X2X3 


X2X3 



TABLE I: The complete and the approximate forms for the functions F'^^ in all channels under 
the normal hierarchy. Note that the approximation for F'^^ in < vl ^ is the same as that 
with Xj + yj ~ and is omitted. Note also that F'^^ = F^°^ . 



where n,f = V^^fi,V*^V;, (Eq. (i)), and 

= ^^JL/AE, (35) 

with L = baseline length. We can rewrite the probability functions in terms of the physical 
observables {x,y). Let us write, for a (3, 

P{iya ^ ly/s) = -4[F2f sin^ $21 + Fgf sin^ $31 + Fgf sin^ $32] 

- 8Jsin$2i sin$3i sin$32. (36) 

Using Eq. (Q, F^'^ can all be expressed as quadratic forms in (x, y). They are listed in Table 

I. The functions F^^^ can be further simplified by using the approximate — t symmetry, 
Xj + ?/j ~ 0, for all A. In addition, with normal hierarchy, xi <C 1 for Ai < A < A^. These 
approximate results are also listed in Table I. Despite the fact that 0:3 ^ 1 for ^ A < Aj, 
terms containing X3 ^ |Vi3p/2 are kept so that the physical potential can be explored. 
Finally, it is noteworthy that the term X2X3, according to Eq. (l23l) . has a resonance behavior 
near A ~ A/^. This is a distinctive feature that can be exploited by proper choices of 
parameters in an experiment. 
For a = /3, we write 

P{i^a ^Va) = l-Y.Pi^^o.^yp) = l + ^Y. sin^ (37) 

a^l3 j>i 

We list F??" in Table II. 

Our results may be compared to formulas in terms of the "standard parametrization" 



381, given, e.g., in Kimura et al. [12|. The relations between {x,y) and the "standard 



parametrization" are given by 

J = Ksiiup 



K = S12C12S13C13S23C23 

222 

Xi = C;^2Ci3C23 — -f^ COS9? 
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complete with Xj + yj Aj < A < Ag 



^21 


-X1X2 + Xiy2 + X2yi 


- yiy2 


— 4xiX2 


< 1 


^31 


-X1X3 + xiys + xsyi 


- ym 


A 

— 4X1X3 


<C 1 


-'^32 


-X2X3 + X2y3 + a:3y2 


- y2y3 


-4X2X3 


-4X2X3 




-X1X3 + j;3y3 + xiy2 


- y2yz 


-X1X2 - (X3/2) 


-X3/2 




-0:2x3 + j;3yi + X2y2 


- yiy2 


-X1X3 - (X2/2) 


-X2/2 




-X1X2 + Xiyi + X2y3 


- yiy3 


-X2X3 - (xi/2) 


-X2X3 




-X2X3 + a;2yi + X3y3 


- yiy3 


-X1X2 - (X3/2) 


-X3/2 




-X1X2 + a;2y2 + a;iy3 


- y2y3 


-X1X3 - (X2/2) 


-X2/2 


J71TT 

^32 


-X1X3 + a;3y2 + xiyi 


- yiy2 


-X2X3 - (xi/2) 


-X2X3 



TABLE II: The complete and approximate forms for the functions F-^^ in all channels under the 
normal hierarchy. 



X2 = 


222 

^12^13^23 COS 


X3 = 


2 22_|_2 2 2_|_ 
'^12'^13'^23 + '^12'^13'^23 + 


yi = 


2 2 2 
-Cl2Ci3-S23 - A COS 


y2 = 


2 2 2 
-^12Cl3C23 - ^ cosy? 


1/3 = 


2 2 2 2 2 2 
'^12'^13"^23 C;^2'^13'^23 



s 



f-Kcos(p 

13 



-.2 



+ l^^f Kcosip (38) 

where Sij = sin 6ij, Cij = cos 6ij, and (p is the Dirac CP phase. It can be shown that the 
functions here in terms of {x,y) are simply ReJ^-'^ in Eqs. (15-23) of Ref. 13], and 



the resultant probabihty functions are identicaL Eq fl551) also offers some insight on the A- 
independence of the approximate fi — t symmetry. It is seen that the conditions Xi + = 
are fulfilled if 1) cL_= s^^, and 2) S12C12S13S23C23 cos = 0. The behaviors of Sij were given 
in Fig. 6 of Ref. [32|- While 5^3 = 1/2 is almost independent of A, S13 = for low A, and 
C12 = for high A. They combine to validate conditions 1) and 2), for all A values. The 
other possibility is that cos ip = 0. Here, (p itself is largely A-independent because of the 
matter invariant sin p sin 2^23 35|| . 

Exact fi — T symmetry was studied earlier by Harrison and Scott [2^. Their formulation 
uses the mixing matrix V (with specific choice of phases), while our results are in terms 
of rephasing invariant (and observable) variables, making it possible to calculate transition 
probabilities directly. In addition, by comparing with the exact formulas in Table I, one can 
quickly compute corrections to the presumed exact symmetry. 



V. APPLICATIONS TO THE LONG BASELINE EXPERIMENTS 

The unique features of the {x,y) parametrization can be used to facilitate, e.g., the 
analyses of the LBL experiments. As an example, let us consider the probability P(t'e — ^ J^^i) 
explicitly. According to Table I, with the approximation Xi + yi = 0, 

P(z/e z/^) = 8[a:;ia;2 sin^ $21 + X1X3 sin^ $31 + a;2a;3 sin^ $32] 

- 8J sin $21 sin $31 sin $32, (39) 
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FIG. 2: For L = 2540 km, the resonant location of 8x2X3 and the peak of the oscihatory factor 
sin^ $32 do not coincide, and the resultant probability P ~ 8x2X3 sin^ $32 is suppressed. Note that 
the probability P ~ 8x2X3 sin^ <I>32 is shown here for a check of the qualitative property at high 
energy. The large, fast oscillating probability near the low energy is not seen here because the X1X2 
and X1X3 terms are ignored in P ~ 8x2X3 sin^ <I>32. 





with J = ±1/2x1X2X3. Using the solutions in Eqs. fl22ll23p . it is straightforward to infer 
the behaviors of P(i^e t'^t)- In the following, let us focus on the region of high A values 
{Ai < A < Ad). Here, xi ^ 1 so that (excluding the case $32 -C 1) 

P{j^e v^) ~ 8x2X3 sin^ $32. (40) 

It is useful to examine the qualitative properties of X2X3 and sin^ $32 separately. If the mass 
hierarchy is normal, the solutions in Eq. ( |23|) suggest a higher resonance for X2X3 at Ah — 
{qh/pl)Ao ^ Ao, where Ao/5o ~ 32. With A/60 ^ [p/(g/cm')][(E/GeV)], So ^ 7.6 x lO^^ 
eV^, and p ~ 3.0 g/cm , the location of resonance A^ corresponds to an energy ~ 10 
GeV, which is independent of the baseline length. Eq. fl40p shows that, in the high A region, 
P(z/e — > u^) (= P(t'e i^t)) is two-flavor like. However, it does not mean that the three- 
flavor problem is reduced to a single two-flavor problem. This is because the probability 
P(z/^ — > z/^), according to Table I, would have contributions from all the $jj's. 

As an illustration, we show 8x2X3, sin^ $23, and P(z/e — )■ f^) ~ 8x2X3 sin^ $23 as functions 
of E in Fig. 2, with L = 2540 km. It is seen that a resonance for 8x2X3 occurs near ii^ = 10 
GeV as expected. However, the smallness of sin^ $32 near ii^ = 10 GeV suppresses the 
probability even if 8x2X3 is at a resonance. On the other hand, the probability at the first 
peak of sin^ <I>32 (near E = 3.5 GeV) also gets suppressed by the smallness of 8x2X3. As 
a result, a significant flavor transition only occurs when L is adjusted so that the peak of 
sin^ $23 is located near the resonance of 8x2X3. 
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FIG. 3: The probability P{i'e ~^ ^n) as a function of E under the normal hierarchy (left) and the 
inverted hierarchy (right). It is seen that P{Li = 7500 km) S> P{L2 = 750 km) near the first 
peak under the normal hierarchy, while P{Li = 7500 km) ^ P{L2 = 750 km) <C 1 for the inverted 
hierarchy. Note that the chosen baseline Li = 7500 km leads to the first peak at £^ = 8 GeV. 



P (Ve^Vj,) 
0.5 



P (Vp^V,, 



Li=7500 km (solid) 
L2=750 km (dashed) 

normal 




E/GeV 



Li=7500 km (solid) 
L2=750 km (dashed) 

inverted 




E/GeV 



The first maximum of sin^ $32 occurs if L/E is properly chosen: 

For the first maximum to coincide with the resonance of .X2.X3, the value of A32 is taken at 
Ah. A32/50 ^ 2^f\V^Ao/6o ^11. It leads to (L/km)/(E/GeV) ~ 10''^ using the current 
upper bound \V13\l ~ 0.03. One concludes that if the mass hierarchy is normal, an extra 
long baseline (L ~ 10"^ km) can lead to a greatly enhanced probability for the neutrino 
beam near ~ 10 GeV, at which energy both 8x2X3 and sin^ $32 reach the maximal values. 
The probability will be suppressed when L starts to vary and sin^ $32 moves away from the 
maximum. Note that for the maxima of a:2a;3 and sin^ $32 to coincide near ~ 10 GeV, 
the baseline L and the undetermined l^islo are related by {L/ km) {\V13\0) ~ 2.54 x 10^. 

On the other hand, since 8x2X3 does not go through the higher resonance under the 
inverted hierarchy, the probability is in general suppressed even if sin^ $32 reaches its max- 
imum. One further concludes that under the inverted hierarchy, the transition probability 
remains small and is insensitive to variation of the baseline length L. 

Thus, if the mass hierarchy is normal, one would expect to observe sizable probability 
difference at high energy for experiments involving two baselines with sizable difference 
in length. On the other hand, the probability would be small and nearly independent 
of the baseline at high energy if the mass hierarchy is inverted. We show in Fig. 3 the 
probability function under both hierarchies for two arbitrarily chosen baselines. Note that 
the peak locations and the peak values vary as L. It is seen that for the normal hierarchy, 
P{Li = 7500 km) > P{L2 = 750 km) near the first peak is expected, while P{Li = 7500 
km) ^ P{L2 = 750 km) ^ 1 if the mass hierarchy is inverted. This result may provide useful 
hints to the determination of the mass hierarchy. Note that the probabilities can be deduced 
if the details of the experiments are considered. If the neutrino energy can be reconstructed 
accurately from the secondary particles involved in an experiment, the observed spectrum 
will tell how the magnitude of the transition probability plays a role. On the other hand, if 
reliable measurement of the energy spectrum is not available, a collection of the event rates 
should also be useful in comparing the probabilities. 

Another possible application is to look for both P{i'e — > u^) and P{Pe ^ix) for a 
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FIG. 4: With the basehne L = 5000 km, P{i^e ^^i) and P{i>e ~^ i^/x) are compared under both the 
normal and the inverted hierarchies. It is seen that near the first peak, P(fe i'^j)/P{pe u^) S> 1 
if the hierarchy is normal, while P{ve ^n)/P{T^e I'n) ^ 1 if the hierarchy is inverted. Note 
that A = 0.02 and the current upper bound for |Fi3|, e = 0.17 have been used. 




single, but very long baseline. Since the P's only go through the higher resonance under the 
inverted hierarchy, one would expect to observe in the vicinity of the peak either P(z/e — )■ 
v^)/P{i'e ^n) ^ 1 if the hierarchy is normal, or P(z/e — )■ z/^)/P(z/e ~^ ^'^) ^ 1 if the 
hierarchy is inverted. We show an example in Fig. 4. Note that although the peak value of 
the probability varies with the baseline length, the relative and qualitative features of the 
above observation remain valid for a chosen baseline. 

VI. CONCLUSIONS 

Neutrino transition probabilities are usually given in terms of the simple expression 
i^aiV/SjVaj^lSi) ^ although the individual V^i^s are not directly observable. When one rewrites 
them using physical observables, such as those in the "standard parametrization" , the re- 
sulting formulas arc often very complicated. It is thus not easy to obtain general properties 
of these probabilities in experimental situations. In this paper we express the probabilities 
as functions of rcphasing invariant parameters. In addition, we incorporate the /i — r sym- 
metry, valid (approximately) for any value of the induced neutrino mass (A). The resulting 
formulas are very simple, and are listed in Tables I and II. They offer a quick quantitative 
assessment for any physical process at arbitrary A values. As an illustration, we analyzed 
the probability P{ue — )■ z/^), with emphasis on its dependence on E, L, and y^jv^Tjl- By 
changing the value of E and L in various LBL experiments, one can hope not only to test 
the theory used to estabhsh P(fQ, — )■ up), but also to help in the efforts to determine the 
unknown parameter |Vi3|q. 
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